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Abstract 


Evolutionary algorithms have been frequently used for dynamic optimization 
problems. With this paper, we contribute to the theoretical understanding of this 
research area. We present the first computational complexity analysis of evolu¬ 
tionary algorithms for a dynamic variant of a classical combinatorial optimization 
problem, namely makespan scheduling. We study the model of a strong adversary 
which is allowed to change one job at regular intervals. Furthermore, we investi¬ 
gate the setting of random changes. Our results show that randomized local search 
and a simple evolutionary algorithm are very effective in dynamically tracking 
changes made to the problem instance. 


Optimization problems in real-world applications often change due to a changing 
environment. Evolutionary algorithms, ant colony optimization and other bio-inspired 
search heuristics have been frequently applied to dynamically changing problems. 

An important approach to gain a theoretical understanding of evolutionary algo¬ 
rithms and other types of bio-inspired computation methods is the computational com¬ 
plexity analysis of these algorithms. During the last 20 years, a large body of results 
and methods has been built up. This includes the development of methods for the anal¬ 
ysis of bio-inspired computing mss and results for some of the best-known 
combinatorial optimization problems such as the traveling salesperson problem 112111 . 

□ll^L and makespan scheduling 122112311 as well as different multi-objective 
olH. These studies often consider the algorithms called Random- 


set cover 
problems 

ized Local Search (RLS) and (1 h- 1) EA, which we also investigate in this paper. Al¬ 
though these algorithms seem to be relatively simple, it should be noted that upper 
bounds on the expected optimization time of these algorithms can often be translated 
to population-based evolutionary algorithms with more complicated variation opera¬ 
tors, e. g., crossover by increasing the upper bounds by only a linear factor with respect 
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to population and problem size 101 ■ We refer the reader to |[ll[Tl , 16 ] for comprehensive 
presentations of this research area. 

In recent years, the computational complexity analysis of these algorithms on d\ 


namically changing problems has gained increasing interest |0, 10 14, 15, T0 2| 


We study one of the classical combinatorial optimization problems, namely makespan 
scheduling on two machines. We consider RLS and (1+1) EA and analyze how they 
are able to keep track of changes that occur to the processing times of the given jobs. In 
our investigations, we examine two models of dynamic changes where in each iteration 
at most the processing time of one job can be changed. In the adversary model, an ad¬ 
versary is able to change the processing time p,- S [E,t/] of an arbitrary job i, possibly 
repeated at regular intervals. First, we show that even for very frequent and arbitrary 
changes, the algorithms are able to obtain solutions of discrepancy at most U frequently 
during the run of the algorithm. Afterwards, we show that RLS and (1+1) EA can main¬ 
tain solutions of discrepancy at most U if the period of changes is not too small. In the 
random model, processing times are from the set {1,..., n} and an adversary is able to 
pick the job i to be changed. The processing time p, of the chosen job is undergoing a 
random change and is either increased or decreased by 1. For the random model, we 
show that the (1+1) EA obtains solutions of discrepancy 0{logn) in time 0{n^\ogn) 
regardless of the initial solution and that the expected ratio between discrepancy and 
makespan is at most 6/n at least once in a phase of 0[r?l^) iterations. 

The outline of the paper is as follows. We introduce the dynamic makespan problem 
and the algorithms under investigation in Section [T] Our analyses for the adversary 
model is presented in Section |2] and the random model is investigated in Section [3 
Finally, we finish with some conclusions. 


1 Preliminaries 

We investigate the performance of randomized local search and a simple evolutionary 
algorithm for a dynamic version of the classical makespan problem. Given n jobs and 
their processing times p, >0, !</<«, the goal is to assign each job to one of two 
machines M\ and M 2 such that the makespan is minimized. A candidate solution is 
given by a vector x G {0,1}", where job / is assigned to machine M\ if x, = 0 and 
assigned to machine M 2 if x, = 1, 1 < / < n. 

The makespan of a candidate solution x is given by 

PA; I 

and the goal is to find a solution x* which minimizes /. We denote by \Mj \ the load 
of machine j = 1,2. We consider the dynamic version of the problem where exactly 
one job changes. We will also allow such changes to be repeated at regular intervals. 
We assume p, G [L,U], 1 < / < n, where L is a lower bound on the processing time of 
any job and U is an upper bound. We denote hy R = U/L the ratio between upper and 
lower bound. 

Randomized local search (RLS) (see Algorithm[Tll starts with a candidate solution x 
and produces in each iteration a new solution y by flipping one randomly chosen bit of 


f{x)=maxl £p,(l--^0:£ 


. '=1 


i=l 
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Algorithm 1: RLS. 


c 

loose X G {0,1}"; 

'hUe stopping criteria not fullfilled do 
y ^.r; 

flip one bit of y chosen uniformly at random; 

_ if f{y) < fix) then .r ^ y 


Algorithm 2: (1+1) EA. 

c 

loose X G {0,1}"; 

'hUe stopping criteria not fullfilled do 
y x\ 

flip each bit of y independently with prob. 1/n; 

_ if fiy) < fix) then .r ^ y 


X. (1+1) EA (see Algorithm|2| works with a more flexible mutation operator which flips 
each bit with probability 1/n. The two introduced algorithms are standard benchmark 
algorithms in the area of runtime analysis of evolutionary computation ililllllQ. 
While evolutionary algorithms usually work with a larger population and potentially 
also a crossover operator, usually positive statements on (1+1) EA transfer to elitist 
population-based evolutionary algorithms by losing only a polynomial factor depen¬ 
dent on the problem and population size |13l- This holds for all results obtained in this 
paper as well as long as there is in each iteration an inverse polynomial probability of 
selecting each individual of the parent population, selection does not accept worsen¬ 
ings of the worst fitness value from the population, and only the variation operator of 
(1+1) EA is applied. 

We study the runtime behaviour of RLS and (1+1) EA on the introduced dynamic 
makespan scheduling problem and their ability to obtain solutions of good discrepancy. 
Eor our theoretical investigations, we do not consider any stopping criteria and measure 
runtime by the number of iterations of the while-loop to achieve a solution of desired 
quality. The expected number of iterations is referred to as the expected time to reach 
the desired goal. In our investigations, we denote by 



the discrepancy of the solution x. We will study the expected time, for different scenar¬ 
ios, until RLS and (1+1) EA have produced solutions of small discrepancy. 

We state an important property on the number of jobs on the fuller machine (i. e., 
the heavier loaded machine, which determines the makespan), which can easily be 
derived by taking into account the upper (U) and lower (L) bound on the processing 
times. 
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• Every solution has at least |'(P/2)/t/] > |’(n/2)(L/t/)] = \{n/2)-R jobs on 
the fuller machine. 


2 Adversary Model 

In this section, we consider the case of a strong adversary. In one change, the adversary 
is allowed to pick one job i to be changed and is able to choose an arbitrary new 
processing time pi G [L,U]. 

2.1 Obtaining a discrepancy of at most U 

We start our analysis by presenting upper bounds for RES and (1+1) EA to obtain a 
discrepancy of at most U from any starting solution. 

2.1.1 RLS 

We first consider RLS and show that the algorithm obtains a solution of discrepancy 
at most U in expected time 6>(nmin{logn,logR}). This bound holds independently 
of the initial solution and the number of changes made by the adversary. The only 
requirement is that the adversary makes at most one change at a time. The proof uses 
the fact that for RLS the number of jobs on the fuller machine does not increase until 
the fuller machine switches. 

Theorem 1 The expected time until RLS has obtained a solution x with d{x) <U is 
(9(nmin{logn,logR}) independently of the initial solution and the number of changes 
made by the adversary. 

Proof We assume that we are starting with an arbitrary solution assigning the jobs to 
the two machines. Let x be the current solution and consider in each point in time the 
fuller machine. The number of jobs on the fuller machine does not increase as this 
would lead to a larger discrepancy. 

We claim that if the fuller machine switched (either by moving a single job or by 
a single change of the adversary) then a solution of discrepancy at most U has been 
obtained in the step before and after the switch. Note that moving one job to another 
machine changes the load on each machine by at most U and that the adversary can 
change the load on each machine by at most U — L. So, the step switching the fuller 
machine (accepted or rejected) reduces the load on the fuller machine from P/2 + a, 
where P = YIi=\Pi^ to P/2 — j3 where a + j5 <U. This implies min{a, j3} <U/2 and 
therefore a discrepancy of at most U directly before and/or after the fuller machine 
has switched. Note, that such a step is only accepted by RLS iff j3 < a and that a 
discrepancy of at most U has been obtained if the step is accepted. On the other hand, 
the case a < j3 which is rejected by RLS implies a discrepancy of at most U before the 
switch. 

The fuller machine has at least [(n/2) • R^^~\ jobs. Let k be the number of jobs 
on the fuller machine. Then the probability to reduce the number of jobs on the fuller 
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machine is ^ and the expected waiting time for such a step is n jk. Summing up, the 
expected time to switch the fuller machine is at most 

A n 

A:=max{[(n/2)-t?“*],l} 

We have two cases. If /? > «/2, the sum is at most «//„ = 6>(nlogn), where H„ is 
the n-th Harmnoic number. If /? < n/2, the sum is at most nlnn + 1 — nln(n/(2/?)) = 
0{n\QgR). Altogether, after at most (9(«min{logn,log/?}) steps a solution of discrep¬ 
ancy at most U has been obtained. □ 

2.1.2 (1-Hl) EA 

In Theorem [T] we exploited that accepted steps of RLS cannot increase the number 
of jobs on the fuller machines. In contrast, the (l-tl) EA may move few big jobs 
from the fuller to the emptier machine and many small jobs the other way round if 
the accumulated effect of the step decreases the discrepancy. Such multiple-bit flips, 
which may increase the number of jobs on the fuller machine, arise in a similar way 
in the analysis of the (l-tl) EA on linear functions, where they complicate the analysis 
considerably lE^I . However, it is also known that the number of incorrectly set bits in 
the (1 h- 1) EA (corresponding to the number of jobs on the fuller machine) has a drift 
towards 0. We are going to show that this drift leads in time to the situation 

that the fuller machine switches, which was analyzed in Theorem[T] We cannot show 
the bound 0{n\ogn) using the advanced drift techniques from ll^ since the dynamics 
of the job sizes do not allow us to use the potential function from the literature. 

Theorem 2 The expected time until the (1+1) EA has obtained a solution x with d(x) < 
U is 0(n^l'^') independently of the initial solution and the number of changes made by 
the adversary. 

Proof We start with a given search point xq, where the time index w. 1. o. g. is 0. 
W. 1. o. g.. Ml is the fuller machine w. r. t. xq. We write 4 to denote the load of Mi 
after t steps. Now, let T denote the first point in time where i, < P/2 + U /2. At this 
time. Ml might still be the fuller machine, which implies a discrepancy at most U. Only 
IHt < P/2 —1//2, the discrepancy is greater than U. Note that Ij-i > P/2 -f t//2 and 
each job size is at most U. Each step resulting mlj <P jl — U /2 must flip at least 2 bits 
and can be converted into a step resulting in fj- G [P/2 —1//2, P/2 -f t//2] by condition¬ 
ing on that a certain subset of bits do not flip. Note that the step defined by the stopping 
time T may be required to flip already more than one bit to reach < P jT + U or 
even no bits may flip at all if the adversary is responsible for reaching Ij <P jT + U jT', 
in the latter case, already discrepancy at most U has been obtained. Note also that 
flipping bits in addition to the ones required to reach Ij < P/2 -f t//2 may result in a 
rejected step. If we condition on the step flipping as few additional bits as possible, we 
are guaranteed to enter the interval [P/2 — t//2,P/2 -f t//2] for the load of the fuller 
machine, resulting in an accepted step. The probability of not flipping a certain subset 
of bits is at least (1 — 1/n)" > Hence, if the step leading to time T flips more than 
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the required bits, we repeat the following analysis and increase the expected time by a 
factor of at most e^. 

We denote by (xf) the number of jobs on M\ with respect to X;, the current search 

point after f steps. Based on this, we define the potential function 


d{xt) : = 



iff < T 
otherwise. 


Hence, the potential function reflects the number of jobs on the fuller machine before 
time T and is set to 0 afterwards. As we have argued, the discrepancy at time T is at 
most U with probabiblity at least 

The aim now is to bound E\T], which is achieved by bounding the expression 
E\d{xt) —d{xt+\) I Xt\t < T] from below and performing drift analysis. In what fol¬ 
lows, we use the notation Xt :=d{xt). 

Since it is necessary to move at least one job from the fuller machine to change 
the (i-value, which happens with probability at least 1 /(en) for each of these jobs, and 
each job on the emptier machine switches machine with probability at most 1 /n, we 
get the bound on the drift 


r , 1 X, / n-X,\ X} 

E[X,-X,+^\Xf,t<T]>^[\ - -]=^, (1) 

which is at least 1 / (en^). Hence, despite the fact that the number of jobs on the fuller 
may increase, its decreases in expectation. Since the maximal t/-value is n, we get 
E\T] = Oin^) by additive drift analysis |01- However, the pessimistic process analyzed 
here has already been more closely investigated in the literature. It has been (apart 
from irrelevant differences in details) modeled by a process called PO-EA by lliol] . 
which was recently revisited by 101. Using this analysis, the bound can be improved to 

In the following, we present a self-contained proof of the 0(r?l^') bound using a 
novel potential function that is easier to handle than the one proposed in the literature. 
Intuitively, our potential function exploits that the process mostly moves due to the 
variance of the one-step change (instead of the very small drift) in the regime X, < %Jn 
whereas it is governed by the actual drift ElX, — Xf+i \ Xt] when Xt is above i/n. 

For X > 0, let the potential function be 


xi\n{^/n) + 2 — ln(x)) 
3>-f 


if X < ^/n, 

otherwise. 


We note that g{x) is monotone increasing and continuous on [0,n]. Moreover, the 
derivative satisfies 


8'{x) := 


dx 


In(Y^) -f 1 — ln(x) 

n 


if X < \/n, 

otherwise. 
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and is non-increasing and continuous as well. Hence, g(x) is a concave function. The 
second derivative equals 

ifx<V«, 

^ dx2 \ — ^ otherwise. 

and satisfies g”{x) < — 1/x forx < n. 

By the mean-value theorem, we get for all x < n and for y > 0 that 

g(x)-g(x-y)>yg'(x) >g(x-fy)-g(x). (2) 

Moreover, by developing Taylor expansions of g{x — y) and g(x-l-y) up to terms of 
fourth order, it is easy to see that 

We are now going to analyze the drift of the process defined by T, := g{X[). To this 
end, it is useful to decompose the drift into a positive and negative part. Define 

Ax:={X,-Xr+i)-l{X,+i<Xr} 


and 

A+:={Xt+i-Xt)-l{Xt+i>Xt} 

and accordingly Ajt and Aji with respect to the T-process. Note that £ [A, — A,+i | A;] = 
E [A-^ I a,] —E [A^ I A,] and accordingly for the drift of the T-process. 

Combining this decomposition with dU, we obtain 


£[T,-T,+i I A,] 

>g'(A,)-£[A^|A,]-g'(A,)£[A+|A,] 

= g'(A,)£[A,-A,+i |A,]. 

If A, > plugging in the expression for g'(Xf) and the bound ([T]i yields 

r ,1 n A^ 1 

£ T,-T,+ i |A,]>^.^ = -, 

A, er en 


which does not depend on A,. 

If A, < i/n, we combine the decomposition with Q and get for some value a{Xt) 

that 


£[T-T,+i |A,] 

> (^fl(A,) + £ [A^ I A,] - a(A,)£ [A+ | A,] 

> +a(A,)£[A,-A,+i I A,] 

■^T 

^ E[Ax l^r] 


7 




since /? I -^jt] ^ ^ sccordin^ to O- Hence, we are left with a bound on 
E [A-^ I A,]. Here we again argue that the number of jobs decreases by 1 if one of 
the Xt jobs from the fuller machine moves and no other jobs moves. Consequently, 
£[A^|^r] >|and 

£[y,-F,+i |A,]> 1 
en 

ifXt < ^/n. Together with the bound derived above, we have E\Yt — Yt+\ \ Aj] > ^ for 
every Xt < n. Now, since Fq < 3n'/^, the additive drift theorem yields E[T]< = 

0{rr’/^) as suggested. □ 


2.2 Recovering a discrepancy of at most U 

We now consider the situation where the algorithms have obtained a solution of dis¬ 
crepancy at most U and the processing time of one arbitrary job is changed afterwards. 
We show an upper bound of (9(min{f?,n}) on the time needed to obtain a discrepancy 
of U after this change. 

Theorem 3 Let x be the current solution that has a discrepancy of at most U before 
changing the processing time of a job on the fuller machine. Then, the expected time of 
RES and (1+1) EA to obtain a discrepancy of at most U is (9(min{/?,n}). 

Proof We use multiplicative drift analysis |I3l to show the 0{n) bound and consider 
drift according to the discrepancy d(x). Let P = Pi and Xt be the random variable 
for d{x) of the search point x at time f > 0. With respect to the parameters from 
the multiplicative drift theorem, we have sq < U + (U — L), imin = U and therefore 
■so/'Smin < 2. W. 1. o. g., let 1,..., / be the jobs on the fuller machine and pi,...pfhe 
their processing times. Furthermore let y{i) be the search point obtained by flipping the 
bit i for As long as the current solution x has discrepancy greater than U, 

each of these single bit flips is accepted. We get 

‘ 4 )' 

> - pA > -{P/2 + d(x)/2)) 

= —{P + d(x)) > — d{x). 
en en 

We set 5 = 1/ (en) and get enln(so/smin) < enln2 = 0(n) as an upper bound. 

It remains to show the 0{R) bound. From the previous calculation, we already have 

E[Xt-Xt+i \ Xt]> —{P + d{x))>P/(en). 
en 

Using additive drift analysis, the expected time to reach a discrepancy of at most U 
when starting with a solution x with d{x) < U + (U — L) is 

(U-L)/(P/{en)) < en{U-L)/{nL)=e{R-\) = 0{R). 


E[Xt-Xt+i\Xt\>-- 

n 


^{d{x) - d{y{i)) 


i=\ 
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Altogether the upper bound is (9(min{/?, n}), which completes the proof. □ 

The previous theorem implies that both algorithms are effectively tracking solu¬ 
tions of discrepancy 0{U) if the time where no changes to the processing times are 
happening is at least c min{/?,n}, where c is an appropriate constant. In particu¬ 
lar, changes happening every c'n iterations where c' is an appropriate constant can be 
tracked effectively regardless of the ratio R = U/L. Furthermore, a small ratio R, e. g. 
a constant, implies that very frequent changes (every c”R iterations, c" an appropriate 
constant) can be tracked by RLS and (l-tl) EA. These statements can be obtained by 
combining drift analysis with an averaging argument over a number of phases. Due to 
space restrictions, this analysis is not spelt out here. 


3 Random Model 


We now consider a model with less adversarial power. Dynamic changes are still possi¬ 
ble, but each change is limited in effect. More precisely, we consider a random model 
as common in the average-case analysis of algorithms Q. For simplicity, we con¬ 
sider the model where all jobs sizes are in {1,...,«}; generalizations to other sets are 
possible. At each point of time, at most one job size can be changed by the adversary. 
The adversary can only choose the job to change, but neither amount or direction of 
change. If a job i is chosen to change, then its processing time changes from its current 
value Pi to one of the two values p, -f 1 and p, — 1, each with probability 1 /2. Two ex¬ 
ceptions are made if p, = n, which results in job size n — 1, and if p, = 1, which results 
in job size 2 afterwards. In other words, the size of each job performs a fair random 
walk on {1,... ,n}, with reflecting barriers. With respect to the initial job sizes, we 
consider both arbitrary (worst-case) initializations and the case that the sizes are drawn 
uniformly at random and independently from Then each initial job size is 

(n -f 1 )/2 in expectation. 

It is useful to denote the random processing time of job i at time t by the random 
variable X,(f). It is well known lll3|] that the random walk described by the process 
Xi{t), f > 0, has a stationary probability distribution given by 


Pr(^/(0 = J) 


1 /(2n — 2) if _/' = 1 or y = n 
1 /(n — 1) otherwise 


Hence, the probability values in the stationary distribution differ from the initial uni¬ 
form distribution by a factor of at most 2. It is also well known in the theory of random 
walks that the so-called mixing time (informally, the time to get sufficiently close to 
the stationary distribution) of the considered random walk is 0{n^) steps. Hence, for 
ij G and for any t > cn^, where t denotes the number of changes to 

job Pi and c is a sufficiently large constant, we have 

n n 

for two constants ci,C 2 > 0. Hereinafter, we asssume this bracketing of 2f,(f) to hold, 
i. e., the mixing time has elapsed for every job. 
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The aim is to analyze the discrepancies obtainable in our model. We summarize in 
the following lemma a useful property of the distribution of the processing times Xiit), 
and drop the time index for convenience. Roughly speaking, it shows that there are no 
big gaps in the set of values that is taken by at least one job. 

Lemma 4 Let (/)(/) := \{Xj \ Xj = iAj G {l,...,n}}|, where i € be the 

frequency of jobs size i. Let G := maxjf | 3i: (/)(/) = 0 (/ + 1) = • • • = (/)(/ + £) =0} 
the maximum gap size, i. e. maximum size of intervals with zero frequency everywhere. 
Then, for some constant c > 0, 

Pr(G >£)< nl-'^^. 

Proof Recall that we assume to be close to the stationary distribution, more precisely 
for each i,jG {1,..., n}, we have Pr(Xy = i) > ci /«. By considering disjoint events. 

Cl 

Pr(Xy G + f}) > — 

for each £ <n. 

Then for each £> 1, we get from the independence of the job sizes that 

Pr{yjG{l,...,n}:Xj^{i,...f + £})< > 

which is at most c^ for some constant C 3 < 1. Hence, by a union bound the probability 
that there is an i such that for all y G { 1 ,... ,n}: Xy ^ {i,... ,/ + f} is at most nc^, which 
equals n2^‘^^ for some constant c > 0 . □ 

Hereinafter, with high probability means probability at least 1 — Ofr’^) for any 
constant c > 0. We prove the following theorem stating that with high probability 
discrepancy 0{\ogn) can be reached in polynomial time. Its proof is inspired by the 
average-case analysis from ll2^ . Note also that the theorem is restricted to the (1H-1) EA 
since its proof analyzes improvements made by swapping two jobs between the fuller 
and emptier machine. 

Theorem 5 Regardless of the initial solution, the following claim holds: with high 
probability the time for the (1+1) EA after a one-time change to obtain a discrepancy 
of at most 0{logn) is (9(n^logn). 

Proof According to LemmalU there is for any constant c > 0 a sufficiently large con¬ 
stant c' > 0 such that there is not gap of size G := c'logn or larger with probabiilty at 
least 1 — n • = 1 — n^‘^. In the following, we assume this maximum gap size to 

hold. 

If the current discrepancy is larger than G, then there must be either at least one 
pair of jobs j,f with j on the fuller machine and / on the emptier machine such that 
Xj! < Xj and Xj — Xji < G, or a job j on the fuller machine of size of at most G. To 
see this, imagine that despite the gap size of at most G, there is no such pair as in the 
first case. Then all jobs of size at least G must be on the fuller machine, resulting in the 
second case. 
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Now, in the first case it is sufficient to swap jobs j and / to decrease the discrep¬ 
ancy by at least 1. In the second case, it is enough to move job j from the fuller to the 
emptier machine to decrease the discrepancy by at least 1. In any case, the probability 
of decreasing the discrepancy is at least (1 — 1/n)"^^^ = Since the maxi¬ 

mum discrepancy is 0{n^), the expected number of decreases is also at most 0{n^). 
Multiplying this with the waiting time for an improvement, we have an expected time 
of By a simple application of Markov’s inequality and repeating phases of 

length cn^ for some constant c, it is easy to see that the time is 0{n^\ogn) with high 
probability. □ 

The previous theorem covers a worst-case initialization with all jobs on one ma¬ 
chine, where the discrepancy can be up to Under random initialization, this is 
unlikely to happen, as the following theorem shows. 

Theorem 6 The expected discrepancy of the random initial solution is ®(n^Jn). Under 
a random initial solution, the time for the (1+1) EA after a one-time change to obtain 
a discrepancy of Oifogn) is (9(n^’^log^«) with high probability. 

Proof We prove that the initial discrepancy is &{n^/n) in expectation and (9(ni/nlogn) 
with high probability. From the last property, the statement on the time to obtain a dis¬ 
crepancy of 6 >(logn) follows with the same ideas as in the proof of Theorem |5] if the 
initial discrepancy is estimated with 0{n^Jn\ogn) instead of 0{n^). 

We are left with the proofs on the initial discrepancy. Let K denote the num¬ 
ber of jobs that are initially put on the first machine. Then E\K] = njl but also 
E[\K-nl2W < •\/Var(/r) = &{y/n), where we used Jensen’s inequality and the fact 
that K ^ Bin(n, I /2). Moreover, by the properties of the binomial distribution we have 
that Pr(^r > n/2 -|-ci/n) =£1(1) for some constant c > 0. Altogether, E[\K — n/2\] = 
0( V^). In other words, there are in expectation 0( Vn) more jobs on one machine than 
on the other. 

Each job size is initially uniformly distributed on n} and has expectation 

(n -f 1)/2. By linearity of expectation, the discrepancy is at least &{^/n) (n -f 1 )/2 = 
&{ny/n). This consideration just subtracts the total load of the machine having the 
minority of the jobs from the total load of the machine having the majority (hereinafter 
called machine “majority”). Should this difference be negative, the discrepancy is still 
positive. However, by approximating the sum of the job sizes on machine “majority” 
by a normal distribution with standard deviation 0 (n^/n), one can also see that the 
discrepancy is 0{n^/n) even if machine “majority” is allowed to have a total load less 
than the other machine. Altogether the expected discrepancy is @{n^). 

The statement that the discrepancy is 0{n^/n\Qgn) with high probability follows 
by using Chernoff bounds on the difference in the number of jobs between the two 
machines (stating that there are at most (9(i/nlogn) more jobs on one machine than the 
other), approximating the tails of the sum of the job sizes on the machines by a normal 
distribution and arguing that a deviation of c-Jn\ogn from the mean has probability 

□ 

Finally, we turn to the case that job sizes change frequently. In the extreme case, at 
every point of time one job size is allowed to increase or decrease by 1. Then it seems 
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hard to obtain a discrepancy of (9(logn) as shown in Theorem|5] However, we can 
apply the results from Section|2] noting that the maximum job size is n at any time. In 
relation to the makespan, the discrepancy, which will also be at most n, is negligible. 

Theorem 7 In the model with random changes, the following holds: the expected 
time until the (1+1) EA (RLS) has obtained a solution with discrepancy at most n is 
0(n^/'^) (respectively 0(nlogn}) independently of the initial solution and the number 
of changes. The expected ratio between discrepancy and makespan is at most 6jn then. 

Proof. Since R = U = nm the notation of Theorem[T]and Theorem|2] we immediately 
obtain the first statement of our theorem. To compute the expected ratio, note that at any 
time the sum of all job sizes has an expected value of «(n+ l)/2 and is at least n^/3 + n 
with probability 1 — using the approximation by Normal distribution. In this 

case, the makespan must be at least n^/6 + n/2, and the ratio is at most nj{n^j6 + 
nj2) < 6/n — 3/n. If the sum of the job sizes is less than n^/3 + n, then the ratio is at 
most n/n since all job sizes are at least one. Altogether, the expected ratio is bounded 
fromaboveby 6/n —3/n + 2^^(") < 6/n. □ 

4 Conclusions 

We have shown that randomized local search and evolutionary algorithms are prov- 
ably successful in tracking solutions of good discrepancy for the dynamic makespan 
scheduling. Investigating the adversary model, we have shown that the algorithms ob¬ 
tain solutions of discrepancy at most U every 6>(nlogn) (for RLS) and every 0(r?l^') 
(for (1 h- 1) EA) iterations even if changes are arbitrary and frequent. Furthermore, such 
a discrepancy is maintained if the period of changes is not too small. For the ran¬ 
dom model, we have shown that discrepancies of 0{\ogn) are obtained and that a ratio 
of at most 6/n between discrepancy and makespan is obtained frequently during the 
optimization process. 
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